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' The possibility of ferromagnetic ordering is revisited in the band model. The coherent potential 

approximation decoupling has been used for the strong on-site Coulomb interaction. The driving 
• forces towards the ferromagnetism are the on-site and inter-site molecular fields coming from differ- 

ent Coulomb interactions. Another driving force is the lowering of the kinetic energy with growing 
magnetic moment coming from the dependence of the hopping integrals on occupation of the neigh- 
boring sites involved in hopping. This effect is described by the hopping interaction, At , and by 
what we call the exchange- hopping interaction, t ex . The exchange- hopping interaction, which is the 
difference in hopping integrals for different occupation of neighboring lattice sites, acts in analogous 
way to the Hund's magnetic exchange interaction. The results are calculated for semi-elliptic den- 
^ , sity of states (DOS) and for the distorted semi-elliptic DOS with the maximum around the Fermi 

energy. They show a natural tendency towards the magnetic ordering at the end of the 3d row for 
C/2 ' the DOS with maximum density around the Fermi energy, when the hopping integrals grow with 

the occupation of the neighboring lattice sites. 
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(S| ; 1. Introduction 

The basic model for magnetic ordering of itinerant electrons in solids is the Hubbard model [1]. The largest 
' interaction in this model is the on-site Coulomb repulsion, U — (i,i\l/r\i,i) , where i is the lattice site index. In the 
mean- field approximation, the Hubbard model leads to the well-known Stoner model for magnetism [2] . The Coulomb 
constant U coming out of the Stoner condition for creating ferromagnetism is large, i.e. of the order of bandwidth. 
On one hand it can be justified by the existing strong Coulomb interaction, but on the other hand for such a strong 
interaction, one can not use mean-field approximation. 

This has prompted attempts to treat the problem within a many body theory. The most significant approach, which 
we will mention in here, is the Hubbard I and conventional CPA [3]. Unfortunately, like many other approaches, they 
also failed to bring any type of ferromagnetic ordering [4]. They did not produce the spin-dependent band shift 
necessary for a ferromagnetic ordering which is why the new versions of the conventional CPA are still being created. 

This is why in our previous paper [5] we investigated the itinerant model for ferromagnetism with both single-site 
and two-site electron correlations. We included also the band degeneration into the model, what has allowed us to 
consider the on-site exchange interactions in the Hamiltonian. The modified Hartree-Fock approximation for the 
two-site interactions was used, which gave us the relative spin band broadening of one spin band with respect to 
the other. This was in addition to the shift in position of majority and minority spin bands. Despite the use of a 
traditional CPA approach the qualitatively new results were obtained which brought the constant of the mean-field 
creating magnetism almost to zero. 

Quite recently some authors (see Hirsch's, Ref. [6]) have pointed out that in the transition to the ordered state 
not only the potential energy is decreased but also the kinetic energy may be lowered. This lowering of the kinetic 
energy is something new in magnetism since the whole previous development has assumed that the kinetic energy 
grows during the transition to magnetic state, as oppose to the decrease in potential energy. The balance of these 
two energies resulted in the Stoner criterion for ferromagnetism. 

The decrease of the kinetic energy, in this new approach, comes from the dependence of hopping integrals on 
occupation of the neighboring lattice site involved in hopping. It is described in our paper by two inter-site interactions, 
which are the hopping interaction, At, and the new interaction, which we call the exchange- hopping interaction, t ex . 
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The hopping interaction is defined as At = to — t\, where to is the hopping integral when no other electron is present 
on sites involved in hopping and t\ is the integral for hopping in the presence of one electron with the opposite 
spin on any of the two sites. The exchange-hopping interaction it is the difference of hopping integrals for different 
occupations of neighboring lattice sites given by Eq. (5) below. For some parameters this interaction is decreasing 
the kinetic energy of the system during transition from paramagnetic to the ordered ferromagnetic state. Both these 
interactions are treated in the linear Hartree-Fock approximation. 

To investigate these effects numerically we will use, like in our previous paper [5], the CPA technique with the 
self-energy describing the on-site Coulomb correlation. This approximation yields the increase in the size (capacity) 
and relative broadening for the majority spin band with respect to the minority spin band. This change of the spin 
band shape is additional to the narrowing of majority spin band with respect to the minority spin band, which is 
coming from the modified inter-site Hartree-Fock approximation. 

The Hartree-Fock approximation is applied to small on-site exchange interaction and small inter-site interactions. As 
a result the on-site small interactions merely contribute to the effective Weiss field, I ex , but the inter-site interactions 
contribute not only to the Weiss field but also to the spin band narrowing different for both spin directions. 

Summarizing, we will have now three driving forces towards the magnetic ordering. 

One is the total exchange interaction, which is the sum of different on-site and inter-site Hartree-Fock interactions. 
This interaction shifts the majority spin band below the minority spin band. Intuitively speaking this interaction is 
coming from lowering the energy of parallel spins with respect to antiparallcl spins [see Eq. (11) below] and it follows 
the idea of Slater [7] and Hund [8] . 

The second driving force is the hopping interaction and exchange-hopping interaction, which are lowering the kinetic 
energy, mentioned above. These interactions in the Hartree-Fock approximation also give rise to the shift of spin bands 
in energy, which is proportional to the magnetization [see Eq. (10)]. 

The third driving force is the change in shape of the spin band under the influence of strong Coulomb on-site 
interaction, U, and the inter-site interactions [9,10] (in here we will use only At and t ex interactions). These effects 
of the change in the band shape, helping ferromagnetism, are mathematically described by the so-called correlation 
factors [4,5,11] Kjj for U interaction and Kf, for At and t ex interactions. 

The hopping interaction and the exchange-hopping interaction, treated in the Hartree-Fock approximation, con- 
tribute also to the superconductivity effect within the BCS formalism (see e.g. Refs. [12-16]). In particular the 
exchange-hopping interaction is the driving force towards the d-wave superconductivity [15] and p-wave superconduc- 
tivity [16]. 

The density of state (DOS), previously semi-elliptic [5], is now assumed to be more realistic; distorted-elliptic, which 
has the maximum around the Fermi energy (see Wahle et. al. [17] with their parameter a =/= and close to minus 
one). 

The paper is organized as follows. In Section 2, we have put forward the model Hamiltonian and developed the 
formalism to treat the on-site and inter-site Coulomb correlation. Numerical examples are presented in Section 3 
based on the semi-elliptic DOS for the weak and strong Coulomb correlation. In addition we analyze the influence 
of the asymmetrical DOS on the magnitude of the on-site exchange interaction. On the base of these results, the 
conclusions regarding the appearance of magnetic ordering with growing occupation of the band are drawn in Section 
4. 



where /io is the chemical potential, cl a {ci a ) creates (destroys) an electron of spin a on the ith lattice site, hi G = cf a Ci a 
is the electron number operator for electrons with spin a on theith lattice site, hi = hi a + hi- a is the operator of the 
total number of electrons on the ith lattice site, h a is the probability of finding the electron with spin a in a given 
band, To is the local energy level, U is the on-site Coulomb repulsion and I ex is the on-site exchange interaction. In 
the Hamiltonian (1) we have three explicit inter-site interactions; J-exchange interaction, J'-pair hopping interaction, 



2. The Model Hamiltonian 



The Hamiltonian for the one-band Hubbard model can be written in the form [1] 



H = (To — Ho) J2hi- Yl tfj ( c ta c jv + h - c ) ~ hx ]C n^flia + hiafli-a 




(1) 
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V-density-density interaction. The spin dependent correlation hopping tfj depends on the occupation of sites i and 
j, and in the operator form it can be expressed as 

tij = t (l - n i _ (T )(l - flj-a) + t 1 [hi- a (l - flj-a) + ftj-„(l - Ui-a)] + tzni-aflj-v, (2) 

where to gives the hopping amplitude for an electron of spin a when both sites i and j are empty. Parameter t\ 
gives the hopping amplitude for an electron of spin a when one of the sites i or j is occupied by an electron with 
opposite spin. Parameter ti gives the hopping amplitude for an electron of spin a when both sites i and j are occupied 
by electrons with opposite spin. Quite recently, several authors suggested that the expected relation to > ti > t 2 , 
may be reversed for large enough inter-atomic distances, t < t\ < t 2 (see [13] and [18]). This concept would fit to 
the results of Gunnarsson and Christensen [19], who for the heavier elements (e.g. 3d or 4f) claim growing hopping 
integrals with increasing occupation. 

For the total on-site exchange interaction I ex one can write on the microscopic level the following expression 

I eX = (d-l)(Jin + JL + Vin), (3) 

where d is the number of sub-bands (orbitals) within the same band, J ini J' in and Vi n are the exchange, pair-hopping, 
and density-density interactions between different orbitals within the same atomic site. In the case of the weak 
correlation this interaction is augmented by the Coulomb repulsion U. 

Including the occupationally dependent hopping given by Eq. (2) into the Hamiltonian (1) we obtain the following 
result 

H = (T - no) Z) "* - I ex E n a nia + n i(T hi-a- + V 

i icr icr <ij> 

+ J E C+Cpda'Cja + J' £ (c+C+C^C^ +h. C ) U\ 
<ij>a,<r' <ij> v ' 

- E [*o - A* (hi- a + hj-cr) + 2t ex n i _ <T rL,-_ CT ] (c+ c jcr + h.c.) 

<ij>a 



where 



At = t -t 1 , t ex = to ^ 2 - ti . (5) 



In this form it is quite visible that At and t ex are also the inter-site interactions. The Hamiltonian (4) will be 
analyzed in two steps. For the kinetic part (the terms with the inter-site interactions At and t ex ) the Hartree-Fock 
approximation will be applied. The other inter-site interactions, V , J , and J', also treated in the Hartree-Fock 
approximation, will be nonzero in the equations of this section to have a full approach, but later on in the numerical 
analysis they will be assumed to be zero as to limit the number of free parameters. The role of the inter-site 
interactions, V, J, and J', was already studied by us before [5]. For the strong Coulomb repulsion U the CPA will be 
used. After setting the energy scale at the atomic level To, and performing the Hartree-Fock approximation we will 
obtain 

H M F = ~ ^ fyf ( C tcr C 3° + h - C ) ~ X] ~ M ^) tl ^ + U Y1 niotli-a, (6) 
<ij><7 ia icr 

where t a e ^ is effective spin-dependent hopping integral given by 

t a eff = t b cr , (7) 
with the parameter b a describing the spin dependent change of the bandwidth 

b a = 1 - 2^n_ CT + 2^ (n 2 _ ff - l\ I a I_„) J —^l a (8) 

to to to to 

The parameter I a = {cl a Cj a ) is the average bond occupation for spin a and the quantity n_ CT above is the probability 
of finding the electron with spin — a in a given band. For the weak correlation one can assume that probability h- a 
is equal to the average number of electrons with spin — a, i.e. fi- a = In the case of strong correlation (U >• D) 
probability of occupation of the band with spin —a, n^ ai will depend on which from the split Hubbard sub-bands we 
are in. For the lower sub-band, e w and n < 1 , we assume that = 1 — n_ CT , but for the upper sub-band, e w U 
and n > 1, we have to assume that h- a = n- a . The modified spin-dependent chemical potential fi is given by 
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jj, = Ho — zVn 



(9) 



and M a is the spin-dependent modified molecular field expressed as 

M a = -I ex n„ + 2zAtI- a - 2zt ex I- a n a - zJn a , (10) 

with z being the number of nearest-neighbors. The expression (10) for the total molecular field shows clearly that 
the difference in hopping integrals t ex plays the same role in creating the exchange field (the third term above) as the 
on-site exchange interaction I ex given by Eq. (11). 

In the simple interpretation of Slater [7], one can understand the total on-site exchange interaction I ex as the 
interaction lowering the energy of each pair of parallel spins with respect to the anti-parallel spins according to the 
equation 

I~ = -I + - ■ (ID 



The similarity of the t ex definition [Eq. (5)] with the intuitive definition of exchange interaction I ex given above, and 
the very same way they contribute to the total exchange field [Eq. (10)] will allow us to call the quantity t ex , the 
exchange- hopping interaction. As it was mentioned above we may have either to > h > t 2 or t < t\ < t 2 , depending 
on the inter-atomic distance. In some cases the hopping integral may grow with increasing electron occupation, 
to < t\ < t 2 (see Ref. [13]), since the overlap between nearest- neighbors atomic orbitals grows as they are expanding. 
We will use both these options below in numerical analyzes of the results. 

For the linear decrease or increase of ti with the growing occupation n, the exchange-hopping interaction from Eq. 
(5) is equal zero. But the dependence t(n) most likely is not linear. In the case of d 2 t/dn 2 > we have t ex < 0, and 
for d 2 t/dn 2 < we have t ex > 0. 

After Fourier transform of the kinetic energy in Hamiltonian (6) we obtain 

hiahi-a, (12) 

kcr io 

where the spin dependent electron dispersion relation is given by 

E%=e k b a -n + M a , (13) 
with £ fc being the initial (without interactions) dispersion energy of the electron 

e k = -tol k , 7fe = E ^ Ri - R ^. (14) 

<i,j> 

In the case of strong on-site correlation U the CPA decoupling is used to analyze Hamiltonian (12), and is described 
by the following equation 

(1 ~ n -° ] l + ~Zk(e) + n -°l-{U-K°)F a {e) = °' (15) 
where E CT is the on-site self-energy and F a (e) is the spin dependent Slater-Koster function given by 

'•<•>- 5? Er^F*- (16> 

This function F a (e) can be expressed by the unperturbed function 

(17> 

k 



with the help of the following relation [5] : 



FAs) = U f s-^+M-V , > (18) 
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which becomes the standard CPA relation when there are no inter-site interactions and in consequence b a = 1. 
For the spin dependent electron DOS one may write the usual expression 

P*(e) = --ImF a (s). (19) 

The spin-dependent average occupation number n a is given by 

+00 

n °= j Pa{e)f{e)de, (20) 

— CO 

where f(e) is the Fermi function 

/ff(£)= 1 + expIe-^-M,)]' P= k^f (21) 

For magnetization per atom in Bohr's magnetons we can write that 

m = n a — n_ CT . (22) 

Differentiating Eqs. (22) and (20), with respect to m and assuming later that m — ► one obtains the criterion for 
the ferromagnetic state in the following form 



l = K + 2 



(dM a 



\ dm 



) j p(s)f(e)exp [/3( £ - n)]de, (23) 



where p(e) is the paramagnetic limit of Pa(s), the correlation factor K describes the role of change in band shape for 
creating magnetization and it is the sum of the on-site and inter-site correlation factors 

K = Ku + K b , (24) 

where 

2 + f° dF a (e) 

Ku - j -^7^ f{£)d£ (25) 

— OO 

and 

+00 

2 f dFj^ dby 

Kb = -n lm J -^d^ md£ - (26) 

—00 

At zero temperature we obtain from Eq. (23) using Eqs. (24)-(26), and (10) the following dependence of the critical 
on-site exchange interaction on carrier concentration 

IZ = 1 ~ 7 ^ (ZJ + lAt + ftex) , (27) 
P(£F) 



lAt 



-4zAt 



dl- 



dm 



and 



2zt ex \ I + n 



dl- 



dm 



(28) 



where Io is the average band occupation in the paramagnetic state. The last term on the right hand side of equation 
(27) is the sum of exchange fields coming from inter-site exchange interaction; zJ, hopping interaction, I At, and the 
exchange-hopping interaction; 7 tcx . The parameter of the bandwidth change b a and the modified molecular field M a 
depend on the quantity I a , the average bond occupation for spin a, which can be expressed as 
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+ oo 



h =< c t c 3° >= ex P[* k ( R i - R j)l J f&Sk* (e) de, (29) 

k -oo 

where 5fc CT (e) is the single-electron spectral density 

SUe) - - > « c+ ;c feCT »= -1 [£ _^_ Re gff} [Im ^ (£)]2 . (30) 
On the other side we can get from the kinetic part of the Hamiltonian with spin a, K 17 , that its average value is 

(K°) = -tl n J2 = -Z% fJ I a , (31) 

<i,3> 

or directly from the definition of the average kinetic energy 

fJ-o 

(K°) = J f{E a )E aPa {e )de , D Q = zt , (32) 

where 

E a (e ) = £o&a, £o = £k- (33) 
Comparing expression (31) and (32) we have 

la = y / -E(j p cj (e ) cfeo, A, = ^o^- (34) 

We can compare now Eqs. (29) and (34). They are equivalent when ReE CT (e) and ImS CT (e) => + (wc remind 
here that E a jD a = Sq/Dq = 7^). This means that I a from Eq. (29) is the generalization of the average bond 
occupation for spin a from Eq. (34) to the case of interaction being described in the single site approximation by the 
self-energy 11(e). 



3. Numerical results and discussion 

The features of this new model will be illustrated by showing the dependence of the critical on-site exchange 
interaction on the carrier concentration I^(n). 

For the unperturbed DOS we will use the asymmetrical function [17] 



Po (e) 



D + ae 



c = 



i + VT 



nD 



(35) 



where D is the unperturbed half-bandwidth and a is the asymmetry parameter. At a = this will become the 
semi-elliptic DOS, and for a = 1 we obtain the DOS for the 3-d fee lattice with t' = to/ '4, where to- nearest- neighbor 
hopping and t' the next-nearest-neighbor hopping [20]. 

The Slater-Koster function corresponding to this DOS can be calculated from the following formula [3] 



F (e) 



(36) 



what will result in the equation 



Fo(e) = 



cirDo 
Do + ae 



e 

A3 



e 



-1 + -- 

a 



- 1 



(37) 



G 



At a = this will give the well-known result for the semi-elliptic DOS 



*b(e) = 




which will be used initially. 

In analyzing our model we will concentrate mainly on how the hopping interaction At and the exchange-hopping 
interaction t ex influence the ferromagnetic state. We will assume all other inter-site interactions; J = J' = V = 0. 
Therefore, the bandwidth factor b a [from Eq. (8)], the chemical potential and the spin-dependent modified molecular 
field will simplify to 

b a = 1 - 2^n_ ff + 2^ (fi 2 _ a - I 2 _ a - 7 CT /_ CT ) , (39) 
to to 

H = Mo - T , (40) 

M a = -I ex n a + 2zAtI_ a - 2zt ex I_ <y fi a . (41) 

The sign and magnitude of interactions At and t ex , according to Eq. (5), depend on the hopping amplitudes t , t\, 
and ti- We assume that t\jto = S and tx/t^ — Si. In general these parameters are different and they both fulfill the 
condition S < 1 and Si < 1 what is equivalent to to > ti > ti (see Ref. 13). However, in his paper Hirsch [13] has 
pointed out that for the hydrogen molecule H2 these integrals depend strongly on the inter-atomic distance and for 
the distance large enough we can even have the reverse relation to < <i < £2 ■ The heavier elements (e.g. 3d or 4f) 
posses larger inter-atomic distances, therefore they may have growing hopping integrals with increasing occupation. 
Gunnarsson and Christensen [19] observed such the dependence for 4f transition elements. In analyzing the influence 
of interactions At and t ex on magnetism we will consider both negative and positive values. Taking into account the 
above defined relations between hopping integrals we can write that At = to(l — S) and t ex = to (1 + SS\ — 2S) /2. 

According to the Eqs. (8) and (10) interactions At and t ex create ferromagnetism by changing the relative widths of 
the spin bands, and by shifting them with respect to each other. The Stoner-Wohlfarth criterion tells us that the large 
DOS on the Fermi level also helps ferromagnetism. In the model with symmetrical semi-elliptic DOS the large value 
of the DOS on the Fermi level can be achieved by decreasing the bandwidth. Another important effect associated 
with changing spin bandwidths is the change of the ratio 6 CT /6_ a . The ratio b a /b- a > 1 helps ferromagnetism for 
concentrations smaller then half-filling, but for concentrations larger then half-filling it is the ratio b a /b-u < 1, which 
helps ferromagnetism [5,9]. The third factor helping ferromagnetism is shifting by interactions the band —a to higher 
energies then the band +o\ 




FIG. 1(a) FIG. 1(b) 
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FIG. 1(c) FIG. 1(d) 

FIG. 1. Critical on-site exchange interaction (FIG. 1(a)), inter-site correlation factor Kb/p(eF) (FIG. 1(b)), band shift factor 
lAt (FIG. 1(c)) and DOS reciprocity 1/p(ef) (FIG. 1(d)) versus carrier concentration, at t ex = and different values of hopping 
interaction; At = 0.4t - solid curve, At = 0.2to - dashed curve, At = — 0.2to - dotted curve and At = Oto - dot-dashed curve. 
All curves in units of Do 

In Figs. 1 we illustrate the role of hopping interaction; At ^ and t ex — 0, but in Figs. 2 we have the reverse 
situation; At = and t ex ^ 0. We will start the analysis with the weak correlation (U <C D). For the weak correlation 
one can assume that probability h- a is equal to the average number of spin —a electrons, i.e. = n- a . For the 
weak correlation the on-site correlation factor Kjj is equal zero. Fig. 1(a) show the dependence for different 

values of hopping interaction At and t ex = 0. Analyzing this dependence one can see that the negative value of 
hopping interaction (S > 1) increases the critical on-site exchange interaction 1^. above the Stoner-Wohlfarth level 
(the dot-dashed curve). Positive hopping interaction At (S < 1) depletes significantly this field especially for electron 
concentrations below the half-filling. According to Eq. (27) hopping interaction modifies the critical on-site exchange 
interaction through the correlation factor Kb, the factor /a*, and the change in the DOS p(e). These factors are 
shown in function of electron concentration for different parameters At = and at t ex ^ in Figs. l(b)-l(d). The 
correlation factor Kb is related to the ratio 6 CT /6_ CT . According to the Eq. (27) to enhance the ferromagnetism we 
need positive Kb, which for n < 1 we obtain at b a /b_ a > 1, and for n > 1 at 6 (T /6_ - < 1. Analyzing K^n), shown 
in Fig. 1(b), we see that Kb is positive for small concentrations and is negative for n > 1. This gives from Eq. 
(8) positive b a /b- a > 1 at all concentrations. Factor I&t shown in Fig. 1(c) gives the spin band shift. For small 
concentrations and At > this factor is positive and the band +er is shifted lower in energy then the band —a what 
helps the ferromagnetism. At n > 1 the factor I/\t is negative, in the result the hopping interaction will oppose there 
the ferromagnetic ordering. The dependence of l/p{e) on concentration shown in Fig. 1(d) has the minimum, which 
shifts towards larger concentrations with growing hopping interaction. Decreasing the bandwidth with growing n (see 
Eq. (8)) what increases the DOS causes this effect. 




FIG. 2(a) 



FIG. 2(b) 
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FIG. 2. Critical on-site exchange interaction ( FIG. 2(a)), inter-site correlation factor Kb/p(sF) ( FIG. 2(b)), band shift 
factor I tcx ( FIG. 2(c)) and DOS reciprocity 1/ p(ef) ( FIG. 2(d)) versus carrier concentration, at At = and different values 
of exchange-hopping interaction; t ex = 0.4to - solid curve, t ex = 0.2to - dashed curve, t ex — ~0.2to - dotted curve and t ex — 
- dot-dashed curve. All curves in units of Do 

Fig. 2(a) presents l^ x (n) for different values of the exchange-hopping interaction t ex and At = 0. Analyzing those 
curves one can see that the exchange-hopping interaction t ex does not decrease much the critical on-site exchange 
interaction but rather shifts the minimums (for positive t ex ) towards larger concentrations (n > 1). Interaction t ex 
modifies value of I£ x through the inter-site correlation factor Kb, which is related with the spin dependent band shift 
factor It Bte and the DOS reciprocity l/p(e) (see Eq. (27)). Their dependence on concentration is shown in Figs. 
2(b)-2(d). Those curves show that for small concentrations the influence of t ex on ferromagnetism is relatively weak. 
At concentrations larger then half-filling and for positive t ex we obtain the positive inter-site correlation factor Kb 
and positive factor 7 tex which enhances ferromagnetism. This enhancement is reduced by the decrease of DOS on 
the Fermi, which at positive interaction t ex comes from the increase in bandwidth with concentration (see Eq. (8)). 
To recapitulate the hopping interaction At enables the ferromagnetism for electron concentrations n < 1, while for 
n > 1 its influence is smaller. The exchange-hopping interaction t ex helps ferromagnetism for large concentrations. 
For n < 1 its effect on ferromagnetism gets weaker. 




FIG. 3 FIG. 4 



FIG. 3. Dependence of the critical on-site exchange interaction (in units of Do) on carrier concentration for the weak Coulomb 
correlation U and different values of the parameter S (Si = S); S = 1 - solid curve, S = 0.6 - dashed curve, S = 0.3 - dotted 
curve and S = - dot-dashed curve. 

FIG. 4. Dependence of the critical on-site exchange interaction (in units of Do) on carrier concentration for the strong 
Coulomb correlation U and different values of the parameter S (Si = S); The curves description is as in FIG. 3 



In Figs. 3 and 4 we show the dependence of the critical on-site exchange interaction versus the carrier concentration 



Iex( n ) f° r the weak and strong Coulomb correlation U. For the curves of I^x( n ) shown in Figs. 3 and 4 we assumed 
that S = Si [18]. In effect we obtained At = t Q (l - S) and t ex = i (l - S) 2 /2. 

The curves presented in Figs. 3 and 4 show that the hopping interaction At together with the exchange-hopping 
interaction t ex decreases the minimum on-site exchange interaction necessary for magnetic ordering. At small 
enough values of the parameter S we obtain the ferromagnetic state for some carrier concentration already at the 
zero values of the on-site exchange interaction I ex . The inter-site correlation factor Kt, depends on parameter of the 
bandwidth change b a , which in turn is a function of the occupation probability of the band with opposite spin, h-a-- 
As it was already mentioned, in the case of the weak correlation (U <JC D) one can assume that this probability, 
is equal to the average number of spin —a electrons, i.e. = n- a . In the case of strong correlation (U ^> D) the 
probability of occupation of the band with spin — <r, h— a , will depend on which from the split Hubbard bands we are 
in. For the lower sub-band (e ~ 0), when concentration n < 1 , we assume that h- a = 1 — n— a , but for the upper 
sub-band (e w U) and n > 1 we have to assume that = n_ CT . For both weak and strong correlation we obtain 
minimum of I%£(ri) for concentrations n < 1. This shows the large influence on ferromagnetism of hopping interaction 
for the hopping integrals fulfilling the relation to > t\ > ti. Therefore, the elements showing ferromagnetism at higher 
concentrations (n > 1) should have large difference between S and Si. 

The large DOS, which favors ferromagnetism according to the Stoner-Wohlfarth criterion, can be achieved by 
decreasing the bandwidth as it was already discussed, but also by considering the asymmetrical DOS. To illustrate 
the results we will use later on the asymmetrical function of Ref. [17] for the unperturbed DOS [see Eq. (35)]. 
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FIG. 5. Dependence of the critical on-site exchange interaction (in units of Do) on carrier concentration for the weak Coulomb 
correlation a and different values of the asymmetry parameter; a = - solid curve, a — —0.3 - dashed curve, a = —0.7 - dotted 
curve and a = —0.97 - dot-dashed curve. 

Fig. 5 presents the dependence Iex( n ) f° r different values of the asymmetry parameter a in the case of the weak 
Coulomb correlation U. For these curves we assumed that At — t ex — which corresponds to the hopping integral 
being independent from the occupation of the sites involved in hopping; to = t\ = <2- Analyzing these results and 
comparing them to the corresponding results obtained for symmetrical DOS one can see that the new DOS causes 
strong asymmetry of the calculated function Ie X (n) (/^(n) ^ i££(2 — nj). For the negative values of the parameter 
a, the minimum of critical on-site exchange interaction will appear at concentrations above the half- filling. At 
positive a this minimum will be shifted to the concentrations smaller then the half-filling. 

The use of the asymmetry parameter increases the DOS at its maximum, what allows for lowering the values of the 
critical on-site exchange interaction as compared to the results for the symmetrical semi-elliptic DOS. 
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FIG. 6. Dependence of the critical on-site exchange interaction (in units of Do) on carrier concentration for the weak Coulomb 
correlation U and the asymmetry parameter a = —0.97. The values of the parameter S (Si = S) are below one corresponding 
to to > t\ > t2\ S = - solid curve, S = 0.3- dashed curve, 5* = 0.6 - dotted curve and S = 1 - dot-dashed curve. 

FIG. 7. Dependence of the critical on-site exchange interaction (in units of Do) on carrier concentration for the weak Coulomb 
correlation U and the asymmetry parameter a = —0.97. The values of the parameter S (Si = S) are above one corresponding 
to to < ti < t2', S = 1.2 - solid curve, S = 1.4 - dashed curve, S = 1.6 - dotted curve and S = 1 - dot-dashed curve. 

Fig. 6 shows the dependence of I^x( n ) f° r a = —0.97 and for different values of the parameter S (S < 1 and we 
assume that S\ = 5). One can see that for small concentrations (n < 1) the decrease of S causes decrease in values of 
the critical on-site exchange interaction Small S corresponds to large hopping and exchange- hopping interactions, 
At and t ex , which in this case help the fcrromagnetism. In a result for small concentrations we have similar behavior 
as in the case of the semi-elliptic DOS (see Fig. 3). The difference is that the small DOS on the Fermi level (for 
a = —0.97) requires for ferromagnetism the nonzero value of the on-site exchange interaction I ex . 

For large concentrations (n > 1 ) and small S (S < 0.5) the increase of S causes the increase of the critical on-site 
exchange interaction I£ x or weakening of ferromagnetism. At larger values of S (0.5 < S < l)the increase of S causes 
decrease of When the parameter S is high enough (S > 1, see Fig. 7) we can get at some concentrations the 
ferromagnetic state without the on-site exchange interaction; = 0. Such a value of S corresponds to the negative 
hopping interaction At and positive exchange-hopping interaction t ex . For these values of interactions, at n > 1, both 
the band shift and the inter-site correlation factor K), are in favor of ferromagnetism. In the case of semi-elliptic DOS 
the ferromagnetic state is not created at these n, At, and t ex , since the DOS is too small on the Fermi level. The 
asymmetric DOS (at a = —0.97) has large values of p(ef), when band is close to full filling, this is why drops even 
to zero. Parameter S > 1 corresponds to to < t\ <t<i- Such a relation according to Hirsch [13], can take place for 
elements with large inter-atomic distances, perhaps elements of 3d and 4f groups. These elements show also strong 
asymmetry in DOS with high density on the Fermi level at the end of the 3d or 4f row. The model presented in here 
could be well fitted to describe ferromagnetism in these elements. 



4. Conclusions 

In this paper we analyzed the influence on ferromagnetism of hopping interaction, At, and exchange-hopping 
interaction, t ex . These interactions are nonzero when the hopping integral depends on the occupations of site involved 
in hopping. 

The physical coupling between these interactions and the magnetic ordering comes from lowering the kinetic energy 
of the system during the transition to magnetic state by these interactions. They depend on the following two 
parameters: S = ti/fo and S\ — t 2 /ti. 

The formalism, which was used in here, was the standard formalism for magnetism with correlation effects. [4,11] 
The correlation effects lead to the change of shape of one spin band with respect to another and also to the bandwidth 
change of one spin band with respect to another. Mathematically they are described by the parameters Kjj and K^. 

The kinetic energy parameters (At and t ex ) influence the ferromagnetism also directly (in the spirit of Weiss theory) 
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by contributing to the band shift factors the following quantities: I At and It ex ■ 

In effect we obtain numerous criteria for ferromagnetism, which show several physical features. 
The hopping interaction At helps ferromagnetism for concentrations n < 1 (see Fig. 1(a)). 

The exchange-hopping interaction t ex helps ferromagnetism, but only weakly, for concentrations n > 1 (see Fig. 
2(a)). 

The on-site Coulomb repulsion U treated in the alloy analogy approximation helps the ferromagnetism mostly in 
the middle of the band (see Fig. 4). 

Fig. 5 shows the critical on-site exchange interaction for ferromagnetism for the asymmetrical DOS with maximum 
density around the Fermi energy. Simple Stoner criterion depicts the strong tendency towards magnetism at the end 
of the band. 

The next two figures display the results of critical on-site exchange interaction for the same DOS, when both kinetic 
energy parameters are present and are the function of S — Si; At = to (1 — S), and t ex — t (l — S) 2 /2. One can see 
that when t 2 < t\ < to (S < 1) the ferromagnetism is enabled at smaller concentrations (see Fig. 6, where the curve 
with S = corresponds to the total exclusion of hopping in the presence of another electron). When t 2 > ii > t 
(S > 1) the ferromagnetism is enhanced quite dramatically at high concentrations (see Fig. 7). This situation can 
correspond to the ferromagnetic pure transition metals. 
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